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ABSTRACT 


A nonparametric identification technique for the identification of close 
coupled dynamic systems with arbitrary memoryless nonlinearities has been 
presented. The method utilizes noisy recorded data (acceleration, velocity 
and displacement) to identify the restoring forces in the system. The masses 
in the system are assumed to be known (or fairly well estimated from the 
design drawings). 

The restoring forces are expanded in a series of orthogonal polynomials and 
the coefficients of these polynomial expansions are obtained by using least 
square fit methods. In most mechanical and structural systems, however, the 
restoring forces can be expressed as a sum of two functions - one related to 
only the displacement of the system, the other related only to its velocity. 
The restoring force in this situation is herein referred to as separable. A 
particularly simple and computationally efficient method has been proposed for 
dealing with separable restoring forces. 

The technique imposes no restrictions on the type of external forcing 
functions used to test the system. The identified results are found to be 
relatively insensitive to measurement noise. An analysis of the effects of 
measurement noise on the quality of the estimates is given. 

Several examples have been provided for the identification of close coupled 
linear and nonlinear dynamic systems. The computations have been shown to be 
relatively quick (when compared say to the Wiener Identification method) and 
the core storage required relatively small, making the method perhaps suitable 
for onboard identification of large space structures. 

vi 



Nonparametric Identification of a Class of Nonlinear 
Close-Coupled Dynamic Systems 

I. INTRODUCTION 

The increased importance given to the accurate prediction of the response of 
structures in various loading environments, has led, in recent years, to a 
growing interest in the improvement of methodologies for proper structural 
modelling. Several investigators have worked in the area of identification of 
structural systems so as to extract, from various types of response data, 
improved characterizations of the systems involved [1-6]. 

The identification problem can be viewed in terms of a class of inputs I, a 
class of models M and an error criterion, E, It usually takes the following 
form.: Given the system response (at various locations in a structure) to the 

class of inputs I, identify a member of the class of models M, which minimizes 
some error criterion E. When sufficient a priori information about the 
mathematical structure of the class M to which a particular physical system 
belongs is available, it is often possible to restrict the identification 
procedure to the determination of the various parameter values, which then 
characterize the system. Such a procedure is referred to as parametric 
identification. On the other hand, as often happens for complex structural 
systems, the a priori information on M may not be sufficient, thus requiring 
the identification procedure to be "expanded" to a search in function space. 
This constitutes nonparametric identification and leads to the "best" 
functional representation of the system. The error criterion E usually takes 
the form of a norm of the difference between the system performance and the 
model prediction. 
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Though a large amount of research has been done in the areas of both 
parametric and nonparametric identification [7-16], present day techniques for 
both are, however, deficient when dealing with large structural systems. 
Parametric methods usually require either the solution of matrix Ricatti 
equations or the employment of nonlinear programming techniques. Often, when 
the number of unknowns in the dynamic system exceeds fifty or so these methods 
in addition to being very computationally expensive may also yield inaccurate 
estimates. Nonparametric methods which employ the Vol terra series or the 
Weiner Kernel approach are also expensive computationally and often do not 
provide adequate characterizations of the types of nonlinearities met with in 
mechanical and structural systems [7-19]. For instance, a "cubic spring" type 
nonlinearity would require the determination of third order kernels whose 
computation in practice becomes prohibitively expensive [20] . In addition, 
the Weiner approach uses white noise inputs. It is often extremely difficult, 
if not impossible, to generate large enough inputs of this nature so as to 
drive large (and often massive) dynamic systems in their nonlinear range of 
response. Applications of such techniques to large, nonlinear, multi-degrees 
of freedom systems are few, if any. 

In this paper we present a relatively simple nonparametric approach for the 
identification of a class of multi-degrees-of-f reedom, close-coupled, 
nonlinear dynamic systems (Figure 1). The method has the advantage of being 
computationally efficient and simple to implement on analog or digital 
machines. Unlike the Weiner Kernel approach, it is not restricted to "white 
noise" type of inputs, and can be used with almost any type of test input. 

The choice of the class of models, M, has been governed by its wide usage in 
problems involving the dynamic response of: (1) full scale building 
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structures, (2) layered soil masses [21], (3) mechanical equipment [22,23], 
and (4) machine components and subsystems in, for instance, the nuclear 
industry [ 24,25,26 ]. 

The technique has been illustrated through application to linear and nonlinear 
systems, and it has been shown that even under noisy measurement conditions, 
the method yields good results. An error analysis is also presented to 
indicate the sensitivity of the method to measurement noise. 

II. SYSTEM MODEL 

The model consists of a lumped mass system with masses M^, £ = 1,..,N which 
are connected to one another by the unknown memoryless nonlinear elements 
K^, £ = 1, 2,..,N as shown. It is assumed that the acceleration of the base 

A/ 

of the structure z(t), and the relative accelerations (with respect to the 
base) X (t), £ = 1, 2,..,N of the various masses are obtained from noisy 

Xj 

measurements. The excitation forces f^(t), £ = 1,2,...N are assumed to be 
available and the masses M^, £ = 1, 2..,N, to be either known or fairly well 
estimated from design drawings. 

Further, for the close coupled system shown (Figure 1) it is reasonable to 
assume that the restoring force depends upon the relative displacement and 
the relative velocity between the masses M^ and • Thus we have 

K^(t) == K^(Y£(t), Y£(t)), = 1,2,. ..,N (1) 
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where, 


^^(t) = x^(t) - 


1,2,. . . ,N 


and 


Xjj^l(t) ^ 0. 


The dot indicates differentiation with respect to time. 


Using noisy measurements of the response x^(t), Z - an estimate 

of the restoring force K is required so that the weighted error norm defined 


£ 


by 


e = 1 |K(y,y) - K(y,y) 


(2) 


is minimized with, G, a suitable weighting matrix and K = {K^}. 

The equations of motion governing the system are then given by: 
M^Xj^ + K^[y^,y^] = - M^z(t) + f^(t) 

^S'^2 ~ ’^I'-^l’^l^ " ~ M2’z(t) + f2(t) 


- ’^£-l[y£-l»y£-l^ = 


( 3 ) 


’Si-l’ ^N-l’^N-l’^N-1^ 
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Equation (3), represents a set of n equations with n unknown restoring forces 
K^, £ == The acceleration of each mass and that of the base, as 

well as the external exciting forces are assumed to be measurable. The 
objective of the report is to identify the unknown restoring forces from a 
measured record of time history of the system. We present below two 
alternative methods for estimating these restoring forces. 


Method I ; 

Adding the top £ equations (£ = 1,2,...,N) together at a time, and 
rearranging, we have the N equations. 


K£[y£,yj^] ~ 1,2, ...,N 


(4) 


where 


w^(t) = {-M^[ (z(t) + x^(t) ] + f^(t)} , £ = 1,2, . . , ,N . 

i=l 

Since w (t) contains quantities which are either known or available from 

X/ 

measurements, the unknown restoring forces estimated. 


(5) 


Method II: 


Equation (3) can be rearranged in the following form: 




(4a) 
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where 


ajj^(t) = -Mj^[z(t) + x(t)] + f^Ct) + 


with 


K 


o 


A 


0 


In (4a), the restoring force ^£^yj^>y£^ depends on the acceleration of and 

that of the base, the external force, and the restoring force on the previous 

mass K (y , y„ ). Noting that K == 0, the restoring force on Mass M,, 
i!,-l 5,-1 5,-1 ® o ’ ^ 1’ 

namely estimated first; the remaining restoring forces can 

be obtained by sequentially using equation (4a) for 5, = 2,...,N. 

III. IDENTIFICATION PROCEDURE 


We particularize the identification approach depending on the characteristics 
of the system under consideration. We consider first, systems with general 
memoryless restoring forces where the restoring forces can be expressed as 
general analytic functions of y and y , and, then, systems with separable 

Aj A/ 

restoring forces where the restoring forces K , 5, = 1,2,..,N can be expressed 

5, 

as a sum of two functions, the first being only dependent on y^ and second 
only on y^. 
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Whereas the first restoring force situation is more general, the 
identification procedure in this case requires a larger computing effort. A 
priori information about the system, especially the foreknowledge that the 
restoring forces are separable, can be used effectively in reducing the 
computing effort. 

Due to the similarity of the basic concepts between Method 1 and II described 
in the previous section, we will devote all our subsequent discussions to 
Method 1 unless otherwise stated. The procedure for Method II is similar with 
slight modifications. 

A. General Memoryless Restoring Force: 

Assuming that the measurements x'^(t), £ = 1,2,.. ..,N and z(t) are 
corrupted by Gaussian white noise, we have the measurements 

x^(t) = x^(t) + and 

( 6 ) 

/N 

'z(t) = z(t) + m(t) . 

Noisy measurements of the various quanitities are denoted by hats above them. 
The noise processes may be assumed to have the following characteristics: 
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E[n^(tj^) nj(t2)] = a» (t^^) 6^(l-j) 6^(t^-t2) , i,j = 1,2, ...,N 


E[n^(t) mCt^) ] =0 i = 1,2,...,N, and 


( 7 ) 




The symbol E[ . ] stands for the expected value, 6^^ stands for the Dirac-delta 

function and 6 for the kroneker delta. 

K. 

The relative displacements x^Ctl and velocities Xjj^(t) may be assumed to be 
either obtained from measurements or from successive integrations of the 
measured acceleration signals. Thus 


x^(t) = x^(t) + Pj^Ct), and 


x^(t) = x^(t) + q^Ct) , 


( 8 ) 


where p ^ t ) and q 1 1 1 , £ = 1,2, 

Xj jo 

equation ("4) then transforms to 


,N, are noise processes. The measurement 


y^] = %(t) ^ w^(t) + v^(t), Jl = 1,2,...,N 

where 


(9) 


£ 

V|^(t) = ^ ^ - lL[m(t) + n^(t)] . 
1=1 


( 10 ) 
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The function ^£[y£> y^] can now be expanded in a double series in terms of two 
sets of functions and Assuming that each set is orthonormal with 

respect to the weighting functions g^^ and g 2 > over suitably defined intervals 
and C 2 , we have 


1,2,. ..,N 


( 11 ) 


j=0 i=0 


. £ 

The coefficients a are to be determined so that the error norm 

ij 

A/ A/ A/ A/ A/ g y ^2 

is minimized, say in the least square sense. This yields the estimate 


IJ 




( 12 ) 


^2 *"1 


where the measurements y. and y are used to replace the exact values y„ and 

■^£ ■^£ £ 


The response quantities that need to be measured for estimating a specific 


K [y , y ], 1 < H < N are then x.(t), i - 1,2,. ..,£+1, and z(t). 

A/ A/ “■ 1 
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B. Separable Memoryless Restoring Force: 


If it is assumed that 


A/ A/ A/ A/ A/ A/ 


(13-A) 


with 


R^[0] = D^[0] = 0, £ = 1,2, ...,N 


(13-B) 


then by (4) we have 


Rn[yp] + D„[yp] = Wp(t), £ = 1,2,. ..,N 

Xj Xj A/ X/ x» 


(14) 


Again expanding Rj^[y^] and D^[y^] in orthonormal sets and we get 




XI 

s=0 


(15) 




q=0 


Estimates of b^ and d^ can be similarly obtained by minimizing G£ in the least 
square sense. 
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In the case where equation (13) is applicable, a simpler alternative approach 


may be followed. As x^(t), £ =1,2, ...,N is measured for te(0,T), the 
quantities y„(t) and y„(t) can be obtained through integration, if x^^Ct) and 

A/ 


£ 


£ 


X (t) are not actually measured. Thus times t 0 ( 0 , T) can be found such that 

£ k,£ 


y^Ctk = 0 ; k = 1,2, ...,k^ , £ = 1,2,. ..,N 


(16) 


For each time t^ ^ which satisfies equation (16) the value of ^) can be 


£ k,£ 


/s 

obtained. As the times t will, because of measurement noise, be slightly 

K. j X/ 

different from those at which y equals zero, '^Qty,(t,)] though close to zero 

36 A/ X< ^ 

may not exactly equal it. In fact if £ = tj^ ^ where tj^ ^ is such 

that yj)(tj^ ^) = 0 then using (13), 


dD. dy„ 
Z ' \ dy„ dt 


k, a 


(17) 




k, £ 


Thus 




(18) 


The coefficients b^, s = 0,1,2,...N , can now be estimated by minimizing 

s R 


2^ 


k=l 


R 






V8<='£<'^k,£>> 


s=0 


; £ = 1,2,. ..,N 


(19) 
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Estimates of then require the solution of the normal equations: 



s=l 




( 20 ) 


where 


I 

k=l 


( 21 ) 


A/ 


k=l 


and the quantities y£(tj^ , k= 1 , 2 . . . ,k ^ have been replaced by their estimates 
y (t ) = y„, . By a proper choice of {(j) } (e.g. Chebyshev polynomials) and by 

A/ K ^ X/ a/K. n 

a proper selection of the points y* (actually achieved in practice by 

a/K. 

interpolation), the matrix S . may often be made a diagonal matrix, so that 

s,j 

''P 9 

h ~ C.T. (22-A) 

J .) .1 

where can be thought of as a normalization constant. 

As k^ becomes large, and the measurement noise tends to zero, arbitrarily 
precise approximations to R£[y£]will be obtained by considering a variety of 
excitation inputs. 
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similarly, the set of points t e(0,T) can be found so that y„(t JO = 0, 

P j A/ iO P ^ 

p = 1,2, ; £= 1,2, ...,N, yielding a normal set of equations similar to 

£ 

(20) and (21), with g replaced by g^, replaced by y^ and the functions 

(y. )} replaced by {(|) (y^)}- Again, by a proper choice of } and a proper 
n X/ n 36 n 

selection of points y , the estimate of d can be expressed in the form 

Pk j 




£ 

d. = E. T 
J J J 


(22-B) 


where E is a normalization constant. 

j 


The method outlined above, is schematically illustrated (for noise free data) 


in Figure 2. We begin with the time histories w^^(t), Y^Ct) and y^(t), £ = 

1,2,...,N. The various times t^ k = l,2,..k^, and p = l,2,..p^, at 

which Y£^t) and y£^t) are respectively zero and determined. The values of 

y’^(t ) and y’^(t .) corresponding to these times are obtained (Figure 2a). 

X/ K. , X/ X/ p ^ X/ 

• ^ 

The corresponding restoring forces R [y*(t „) ] and D [y (t )] are found as 

X< K y X/ X/ P y Xj 

w^[tk and w^[t^ respectively. The values of and are then plotted 
versus y^ and y^ respectively (Figure 2b). 


£ 


We mentioned in passing that often sufficient a priori information on the 
nature of some of the restoring forces may be available, i.e. some of them may 
be known to be of the separable type. Noting that for purposes of identifi- 
cation, each is uncoupled from the others (equation 4), a combination of 
the methods presented above can be used - the general case for all the K^'s 
which have general restoring forces and the separable case for all the E^'s 
which have separable restoring forces. 
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c. 


Forced Vibration Testing of Systems 


In the absence of a base motion z(t), the identification of K£ly£>y£J 
can be performed without the need of explicitly knowing (or 
measuring) the forces f^(t) if we specify that 


f±(t) 


arbitrary 


1 < i < SI, and 


i > i 


(23) 


For noise free data, various arbitrary functions f (t), i-£, can be used so 

1 

that arbitrarily accurate approximations of can be found. 


IV. ERROR ANALYSIS 


• ^ 9 ^ Q • • • • 

The estimates a , b^ and d^ obtained by the simple regression analysis 

Ij 1 1 

technique outlined in the previous sections differ from the exact values 
primarily because of the presence of measurement noise. The influence of 

noise may be thought of as affecting; (1) the measured value of w (t) and (2) 

SL 

the estimates of y (t) and y (t). 

i i 

To acquire an appreciation of the effect of measurement noise on the 
estimates, we shall consider here the case where the restoring force is 
separable. Error analysis of equation (12) for the general restoring force 
case, besides being more complex, will not, it appears, yield any additional 
physical insight into the effect of measurement noise on the estimates. 
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Let the discrete time points t „ be utilized where the times t . are chosen 

p,^ p,£ 

/S 

so that, say, y„ft .) = 0, p= l,2...,p, I = 1,2...,N. 

3o p j X/ A/ 


If we assume that the noise in measuring, x (t ) and Z (t ), 

p j Ai p j 36 

p = 1,2, ...,p ; £ = 1,2,, ...,N, has zero mean, is uncorrelated and has a 


constant variance, then by equation (10), 


E[v^(t^ ^) ] — 0, 'V t £ (0,T), l,2,,,.,Nj 


and 


(24) 






*^p,£’ ^q,£ ^ 1,2, ...N, 


where 


a„ = 


2 

O’* ( t n ) 
X p,X 


E 

i=l 


k,2 1,2, 

M. I + 0,. (t J 
If z p,X 


£ \2 

y«. 

1 

1=1 


(25) 


Furthermore, if Xj^(t) and x^it),?- =1,2,...,N+1 are measured, then x^e have the 
following relations: 


and 




with p (t)=q (t)=0. 

N+1 N+1 
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The random variables and 3j^ are assumed uncorrelated, such that for any two 

times t and t .e (0,T), 

p,£ q,£ 


E[a^ (tp^^)] = E[3j^ ^ = 1.2,... ,N, 


^ \ \ (p-q). ^ = 1,2,...,N-1, 




(26-B) 


^^^£ ^^p,£^®£ ^*^q £^^ "" ^ ^ 1,2,...,N-1, and 


The variances of the random errors in measuring (t and (t ^), 

2 2 .* * 

£ = 1,2,...,N are assumed to be cr and o. respectively. The measurements of 

^ , A 2 

z (t 5 ) and z (t „) are also corrupted with noise whose variances are cr and 

P J A/ P j X/ " 

2 . . 

G . . The various variances could be functions of time, 
z 


Let tp ^ be such that y^ (tp = 0, P = 1,2, ...p^. 


(27) 


Define t „ by the relation 

p , X 


*^p,£ ^p,£ ^p,£ 


( 28 ) 


where T n is the error in finding t n. 

p , X, P » ^ 


We 


next expand y^ (tp ^) in a Taylor series giving, 
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tit 


1 2 ,2 

‘‘"<1 


£ 


Pp-^ 


2 , 2 
dt „ 

p,^ 


+ higher order terras . 


(29) 


Using equation (24) and noting relations (26-A) and (23) we get 


✓\ /s 


dye 


^ ~ ^p,5- dt 


P.2. 




(30) 


where the higher order terras in n have been neglected. Thus if 

U «k X/ 






(31) 


For the oscillating system considered, it will generally be unlikely that 

• 

y^(tpj^) and Y£(tp ^) b® both zero, except when the oscillator is executing 

very small amplitude motions, preparatory to coming to rest. If, however, 

• 

y„(t „) = 0, then the next higher term in (29) can be used to estimate T 
A/ p j ^ p j X/ 


Then w (t .) can be expanded as 

p , X. 


^0 (t„ (b n) + a nOtp (t (j) ~ + 


/N 


'"p,r ""ii ’ °p,rn ^"p.A^ dt 


(32) 


P.2 


dt 


+ higher order terms + v„(t ) 

£ p , £ 


and y ( t ) as 
p,£ 
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,H “r 




A 

+ higher order terms + (t^ 


(33) 


Neglecting the higher order terms, equation (33) gives, 


A 



A 




(t J 
p,il 


+ Yc 


A 



(34) 


where , 




.l^SL 










+ 3 


SL 


A 



(35) 


Using equation (22-B), we have 


d. 

3 


E. 

3 


P£ 


E 

•r»^ 1 


E. 

J 


A 





A 



A 



(36) 


where 4^, = gy '<P . > 
J 3 


Taking the expected value on both sides of equation (36) we get. 


''0 

E[dp 


p=i 


A A 




( 37 ) 
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The function ip. can be expanded about yn(t „) to yield, after some 

j X- p , X. 

manipulation, 


E[dJ] 


I 




p=l 


p=l 




+ higher order moment terms 


where 






2 / ='l"2 dT! + y*"e 7,72 + 


'i ^ Ay 


% 


=i 


,27 

d V, 


w 


£ ,'2 

J 




and. 


dip, 

9 

dy. 


'P»^ 


p - i i 2 
p,£ 4 j ^p,£ ^^£ ^- 


3%.: 


3y 


Z t 


p,£ 


(38) 


(39) 


(40) 
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similarly it can be shown that 


I 

E[t>tl CjW„(t,_ „) ())_, (y„(t,^ „)) 


) . 

J 


k,r 


k=l 


^ E[aJ(t,, „)] A,* „ + E[3^t, „)]B,’ 


k,V‘ ^k,£ ■'"k,£ 


k=l 


2 2 

+ E[a^(t^ ^) ] E[g^(tj^ ^^^^k a ^ ^ higher order moment terms. 


with 


A' 


k,H 


w 


d^$. 

3 

£ , 2 


dyr 


-^t 


k,£ 


bV 


k,£ 


d^. ^2 d(|). ^ 

k,£ 2 ( ^ + ^£ ^^£ ^ % ^3 + ^ 


k,£ 


and, 


";.£ 


2 “^ 

1 I 2 .. ^ 

T ^b, „w„ 


4 rk,£™A 3y2 


k,£ 


where 


\,£ : ^ = ®1 '^3 ' 


( 41 ) 


( 42 ) 


( 43 ) 
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Had the signals and been obtained by integrating Xj^^, the errors and 
3^ would have been correlated at various times leading to additional terms in 
equation (40) involving the expected value of the products of (t^ and 


Equations (38) to (43) indicate the effect of the measurement noise on the 

a £ 

expected value of the coefficients b . and d .. We observe that the 

3 3 

estimates are biased, the bias being independent of the noise in the 

measurement of x^(t) and z ^(t) , The bias is however dependent on the noise in 

the measurement of x^^(t), x^(t), ^nd x^^^(t). If the noise in these 

measurements goes to zero, (i.e., ot^'s and s equal zero) all the terms 

except the first on the right hand side of equations (38) and (41) go to zero 

so that the estimates become unbiased. Furthermore, if say, we use the 

orthonormal polynomials {m } or {(b } , greater biases would, in general, be 

3 3 

• • ^ P • • • • • • 

obtained for the coefficient estimates d and b with increasing j. This is 

j j 

because the higher order polynomials oscillate more rapidly thus leading to 
larger values of the dij;./dy„, d^ijj./dy?, etc., which in turn, by equations (39) 
and (41), increase the bias. 


Further, to illustrate the effect of noise in the measurement of the 
acceleration terms, let us assume that ap(t), B/t), £= 1,2,...,N, are zero. 

Xr A/ 

We have then 

^0 

Var [d . ] = E 

;i 


p=l 


(44) 
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which for uncorrelated noise gives 


Pjl 

Var [d*l - ^ E[vJ ' 

P = 1 


( 45 ) 


But we have by equation (25) 



which indicates that the variance of the estimate increases with Z as well as 
with increasing magnitudes of the masses M^, i = 1,2,...,£. Since > 0, ¥Z, 
it follows that 



Thus from equation (43), the variance is more sensitive to noise in the 
measurement of the base motion z(t), than it is to noise in x^(t), Z = 
1,2, ••• ,]SI« 


V. APPLICATIONS 

In this section, a few select applications of the identification technique 
discussed earlier are presented. The responses of the systems considered are 
simulated by integrating the dynamic equations. Noisy measurement records are 
obtained by adding zero mean Gaussian white noise (ZMGWN) to the integrated 
results. 
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Both the separable and the general restoring force situations have been 
illustrated. Method 1 and II (Section II) are applicable to both categories. 
For the separable restoring force case, the linear system (described below) is 
identified by Method I using calculations in single precision. On the other 
hand, the nonlinear systems have been identified by both methods in both 
single precision and double precision. A comparison of the results of these 
methods is reported. The general restoring force case is investigated using 
only computations in single precision. 

Motivated by the simplicity of the method, it was attempted to investigate 
its worthwhileness in a simulated real time environment using a small computer 
with a maximum core storage of 64KB for computations in single precision. For 
double precision work, a bigger core was used. A sector of forty seconds of 
data in each case was analyzed. The digitization rate for the data was taken 
to be 0.04 sec, a rate which would allow the multiplexing of several channels 
using standardly available analog-to-digital convectors. The model used is a 
four degree-o£-f reedom system as shown in Fig. 3. 

A. Separable Restoring Force Case. 

Three different systems are considered in this category, a linear system 
and two nonlinear systems. 

a) Linear System 
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Consider the system shown in figure 3 with the restoring forces given by 


Kj^[y,y] = ^ "" 1,2, 3, 4. 


If the system is linear, then 


R^[y] = 

and, 

Djtyl = djy 


£ = 1 , 2 , 3 , 4 . 


(47) 


The various parameters of the system are shown in Table 1. The system is 
subjected to the swept-sine wave test excitation, 


f^(t) = aj^Sln[u)(t) ] t i = 1,2, 3, 4 


(48) 


where the time dependent frequency w(t) changes linearly on the time interval, 
(0,T) according to the relation. 


a)(t) = + a^t 


(49) 


where. 


n^a^ 

“2 ■ 


(50) 


n^ and ^2 are scaling constants, and Tq is a normalizing time constant. 
Figure 4 shows a segment of the excitation signal (described in Table 1) at 
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each of the four mass levels for a = 2 rad/sec, n = 5, n =40, 


OJ ^ (271/1 ) = , / b^/m = 10, and T = 40 secs. The time scale is shown 

normalized with respect to T^, A short portion of the system response to 
this excitation is indicated in Figure 5. Forty seconds of data (approxi- 
mately 15 times the fundamental period of the system) are used for the 
identification. By digitizing this data at equispaced time increments At = 
0.04, the x^(t), digitized time histories and x^(t) are obtained. 

To study the effect of measurement noise on the identification results, these 
digitized results are corrupted by the addition of zero mean uncorrelated 
gaussian noise. The same noise-to-signal ratio (N/S) is used for each of the 

measurements X , x , x , % - 1,2, 3, 4. The identification results are obtained 

III 

for three different values of the N/S ratio, namely, 0.001, O.Ol and 0.02. 
Whereas the first number represents data of exceptionally good quality, the 
second typically represents the situation pertinent to data available from 
accelerographs, and, the third to what may be referred to as "poor" quality 
data. 


From these 'noisy' measurements, the corresponding 
and w^(t) are calculated for t = iAt, i=0,l,2,.... 


time histories y (t), 

I 


1000 . 


y (t) 

I 
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The functions and are expanded in a series of Chebychev polynomials 
so that , 


N. 


R 




s=0 


s=0 


and 


( 51 ) 






q=0 q=0 


q^£ 


The values of and are chosen to be 4 and 3 respectively. The 

coefficient estimates b and d are obtained (by performing a least square 

fit) by solving the normal equations (Equations (20) and (21)) [27]. To 

improve the quality of the fit [27], the y(tj^ ^ and y(tp £) arrays are 

normalized so that they lie in the interval (-1,1). Using the weighting 

2 — 1/2 Z ^ Z 

functions gi(n) = go(n) = (1-n ) , the coefficients b and d are found. 

For ease of comparison with the exact and D^'s, these coefficients are 


Z 


^ Z 

converted to b and d corresponding to the polynomial expansions (Equation 
s q 


51). 


Figure 6 shows the results of the identification giving the estimates of the 

intermass stiffness (R^) and the intermass damping (D£) as functions of 

relative displacement and velocity respectively. The least square polynomial 

fits are calculated at the various points y(t ) and y(t ) for various 

k,£ k,£ 

noise-to-signal ratios. The exact stiffness and damping are also plotted at 
the same values of y and y for comparison. As seen from the figure, the 
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estimates gradually worsen with increasing values of N/S. The estimated 
coefficients of the polynomials are shown in Table 2(a) and (b) for each R 

and I) , a == 1,2, 3, 4. We observe that, in each case, the estimated 

SL 

coefficients for all except the linear term are small. 


£ 


A measure of the accuracy of the identified stiffness and damping can be 
obtained by defining the root mean square errors (rms) as 


/*-^£ ^£^ 

/ 


, 1/2 




/ 


-. 1/2 



(52) 


where the integrations are carried out over the complete response range of y 


£ 


and y respectively. The rms errors are indicated for each R„ and D and each 

£ 36 36 

N/S ratio in Tables 2(a) and 2(b). 


It is interesting to note that the rms error does not change substantially 
when the N/S ratio changes from 0.001 to 0.01. This is because of the fact 
that for such low values of the N/S ratio the digitization process as well as 
the single precision accuracy of the computations (which leads to round off) 
actually dominates the accuracy of the results. We note from the tables that, 
in accordance with our discussion in Section IV, the rms error increases with 
increasing 'i' values. 
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A, 


A comparison of the predicted response using and and the exact response 
for an excitation different from the test excitation, and comprising a base 
acceleration, z(t), is indicated in Figure 7(a). This base acceleration is 
actually a sample of zero mean Gaussian White Noise (ZMGWN) with a standard 
deviation (a) of unity. The stiffness and damping estimates corresponding to 
the N/S = 0.02 case are used. We observe that the predicted responses, using 
the identification results obtained even under very noisy test conditions 
(N/S = 0.02), and the exact responses are reasonably close to each other. 


The solid lines in Figure 7(b) show the response of the system when mass is 
subjected to an impulsive (delta-function) force of ten units. The predicted 
response of the system, using the identification results obtained for N/S = 
0.02 (Tables 2a and 2b), is indicated by the dashed lines. Again, the 
predicted response matches well with the exact response. 

b) Nonlinear Systems 


Two nonlinear systems have been considered. They represent nonlinearities 
which are often encountered in structural and mechanical systems. The first 
system has nonlinear stiffness and linear damping of the form. 


and 


D r 1 - a. i,^y3 

^1 £ ^3 £ 


°£t^£^ = V£ 


£ = 1 , 2 , 3, 4 


The system description is given in Table 3. This system is identified by 
method I using single precision computations. We note that whereas R , R , R 


28 



represent 'hardening' nonlinearities, R ^ represents a 'softening' 
nonlinearity. The test signal used is identical to that used for the linear 
system described in Table 1. Using = 4 and = 3, b^ and d^ were 
obtained. The estimated functions and are shown, as before, in Figure 
8. Tables 4a and 4b give the estimates for the coefficients of the polynomial 
series representation of Rj^ and D^. The coefficients are obtained via the 
Chebychev polynomial expansion as mentioned earlier. The rms values for 
different N/S ratios are also indicated. It is seen that the identification 
procedure leads to fairly good estimates even when using noisy (N/S = 1/50) 
test data. 


Figure 9(a) shows a comparison between the predicted response of the system 
(using the identification results of Tables 4(a) and 4(b)) and the exact 
response of the system when the system is subjected to twice the amplitude of 
the ZMGWN base acceleration used before (Figure 7(a)). Identification results 
corresponding to the N/S ratio of 0.02 were used. Figure 9(b) shows the 
predicted and actual system response to an impulsive force of ten units 
applied to mass M^. 


Secondly, a system with nonlinear stiffness and nonlinear damping is chosen 
and identified by both method I and method II. The system used is the one 
with 


and 






r’ 1 


Z = 1 , 2 , 3, 4 
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Table 5 shows the actual parameters of the system. Identification of the 

coefficients b and d is done using N_ = = 3 in both methods with the 

~s ~q K u 

test signal defined in Table 1. Since the records of the system are obtained 
through the integration of the dynamic equation of the system. The accuracy 
of computation may be one of the factors which affect the estimated 
functions. An additional double precision calculation, besides a usual single 
precision calculation, is used in both methods. The estimated functions R 


Z 


and done by method I in single precision calculation for different N/S 


ratios as well as the exact functions are plotted in Figure 10. The response 
of the actual system and the identified model done by method I and single 
precision calculation (using results of N/S = .02) subjected to the base 
acceleration of Figure 7(a) and the same impulsive loading used before are 
reported in Figure 11(a) and 11(b). Figures (12-14) shows the estimated 
functions and done by method I and method II in single precision and 
double precision calculations with various N/S ratios as well as the exact 
functions. The figures indicate that the accuracy of the functions estimated 
by both methods is essentially the same. The double precision calculations 
while requiring more computational effort and core space yield marginal 
improvements in the estimates. The results of the estimated functions with 
the RMS errors are shown in Tables 6 (a) and 6 (b). The RMS errors for method 
II are seen to be slightly higher than those for method I. 


B. The Generalized Restoring Force Case. 


Expressing the restoring forces in terms of the orthogonal Chebychev 
polynominals , we have 
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n m 


i=l j=l 


The coefficients a^j are obtained as shown in equation (12) by minimizing 
the error norm in the least square approach. 

Two systems, both nonlinear, have been considered. To compare the general 
restoring force approach with the separable case, the nonlinear system of 
Table 3 is identified assuming no a priori knowledge regarding the 
separability of the restoring force. 

A Chebychev polynomial expansion using 64 coefficients (n = m = 8) is 

• • 

employed. The variables y^ and are normalized to lie between -1 and +1, 
and 600 data points are used for the least square fit. The identified 
coefficients are then converted to monomial bases for ease of comparison with 
Table 3. Table 7 shows the identification results for different N/S ratios 
and the RMS errors involved. It is observed that the identified stiffness and 
damping estimates obtained by this method are in general superior to those 
obtained using the separable restoring force method. This is attributed not 
only to the increased number of data points that are used to perform the fit 
here, but also to the inaccuracies in the separable case that accompany the 
estimation of the times at which the velocities and displacements become zero. 
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The second system considered is described by the relation 




£ = 1,2,. ..,N 


The coefficients a^, b^, c^ , d,^, e^ , f^, £ = 1,..,4, are shown together with 
the identified results for various N/S ratios in Table 8. The identification 
is done using 600 data points. The RMS errors are also presented. Perhaps a 
better way of comparing the identification results with the exact system is to 

A 

compare K [y , y ] and K [y , y„ ] • This is done in Figures 15-17 for various 

36 A/ 36 36 36 

N/S ratios. It should be noted that even though some of the identified 
coefficients differ substantially from those of the actual system, in the 
regime of response considered, the RMS errors are small. It is these RMS 
errors which should be, in reality, considered when assessing the quality of 
the identified results. 


VI. CONCLUSIONS AND DISCUSSION 

A relatively simple nonparametric method for the identification of a class of 
close-coupled nonlinear multi-degree-of-freedom systems has been developed. 

The class of systems is one which is often encountered in the fields of 
mechanical and structural engineering. Identification of arbitrary memoryless 
nonlinearities is possible through knowledge of the accelerations, velocities 
and displacements of the various masses. These quantities are then used to 
obtain by regression techniques the surfaces of the restoring forces as 
functions of the intermass displacements and velocities. 
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A particularly simple and computationally efficient technique is 
illustrated when the restoring force is linearly separable into two functions, 
one of inteirmass velocity and the other of intermass displacement. 

An assessment has been made of the effect of measurement noise on the 
estimates of the coefficients that are obtained from the regression analysis. 
It is found that whereas the biases in the estimated coefficients are 
primarily dependent on the noise in the displacement and velocity 
measurements, their variances are controlled to a good extent by noise in the 
acceleration measurements. 

All the computations related to single precision calculations have been 
performed on a small 16 bit minicomputer, with a 64KB memory. Even under very 
noise measurement conditions, (N/S ratio of 1/50) with only a few terms in the 
seri( 2 s expansion, the identification results yield low rms errors. The 
capability of predicting the response of the system to excitations other than 
the test excitation, by using the results from identification, has been 
illustrated,, As has been observed in other studies [7] accurate estimates of 
the damping are in general more difficult to obtain than estimates of the 
stiffness. Double precision computations while significantly increasing the 
computational effort and the core required (beyond 64KB) did not yield 
substantial improvements in the estimates. 

A drawback of the method is that it can only be used for identifying 
memoryless intermass nonlinear restoring forces. This is so because 
expansions of the type given by equations (11) and (15), where the y„ ' s and 

X/ 
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Yj^' s are treated as independent variables, are only valid if the restoring 
forces are single valued functions of the independent variables. Thus, for 
example, in a bilinear hysteretic system in which the restoring force is a 
multivalued function of and y^, the technique would fail. Alternatively 
speaking, for such systems, one could find a class of inputs which would yield 
incorrect identification. A simple example of such a class of inputs, for the 
bilinear hysteretic case, is the class of inpulsive excitations which cause 
permanent displacements of the system. 

The main advantages of the method are: 

(1) The method is applicable to general meraoryless intermass nonlinear 
restoring forces. 

(2) There is no limitation on the nature of the test excitation that can 
be used for the identification. This is a major advantage over some 
of the other non-parametric methods [20, 21] which often require 
Gaussian White Noise (GWN) excitations. Such GWN excitations are 
difficult to produce in high enough magnitudes in order to drive 
multi-degree-of-f reedom systems, which are often large, in their 
nonlinear ranges of response. 

(3) The computational requirements, both in terms of CPU time as well as 
storage, are very small in comparison with the Weiner Method making 
the method attractive for real time identification [21]. 
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(4) The duration of time over which the data is required to be taken is 
comparatively small compared to other nonparametric techniques [21]. 

(5) The identification results obtained are relatively insensitive to 
measurement noise. The rms errors in the determination of the 
restoring forces increase in general as we move towards the point of 
fixity of the system. 



VII REFERENCES 


1. Eykhoff, P. "Process Parameter and State Estimation", Survey Paper, Proc. 
I.F.A.C. 3rd Congress, London 1966. 

2. Berman, A., "System Identification of a Complex Structure", AIAA 16th 
Structures, Structural Dynamics and Materials Conference, Denver, 

Colorado. 

3. Dale, D.B. and Cohen, R. , "Multiparameter Identification in Linear 
Continuous Vibratory Systems", ASME Journal of Dynamic Systems, 
Measurements and Control, Vol 93, March 1971. 

4. Eshleman, R.A., "System Modeling", The Shock and Vibration Digest, Vol. 4, 
May 1972, p. 1. 

5. Hart, G.C., ed.. Dynamic Response of Structures; Instrumentation, Testing 
Methods and System Identification, ASCE/EMD Specialty Conference, UCLA, 
1976. 

6. Collins, J.D., et al.. Statistical Identification of Structures", AIAA 
Journal, Vol. 12, 1970, pp. 185-190. 

7. Udwadia, F.E. and Jerath, N., "Time Variations of Structural Properties 
During Strong Ground Shaking", ASCE, EMI, Vol. 106, 1980, pp. 111-121. 

8. Beliveau, J.G., "Identification of Viscious Damping in Structures from 
Model Information", ASME, Journal of Applied Mechanics, Vol. 43, June 
1976, pp. 335-339. 

9. Distefano, N., and Todeschini, R. , "Modeling, Identification and 
Prediction of a Class of Nonlinear Viscoelastic Materials", International 
Journal of Solids and Structures Part I, Vol. 9, 1976, pp. 805-818. 

10. Udwadia, F.E. and Shah, P.C., "Identification of Structures Through 
Records Obtained During Strong Earthquake Ground Motion", ASME, Journal of 
Engineering for Industry, Vol. 98, No. 4, 1976, pp. 1347-1362. 


36 



11. Udwadia, F.E. Sharma, D.K. and Shah, P.C., "Uniqueness of Damping and 
Stiffness Distribution in the Identification of Soil and Structural 
Systems", ASME Journal of Applied Mechanics, Vol. 45, March 1978, pp. 
181-187. 

12. Kaya, I. and McNiven, H., "Investigvation of the Elastic Characteristics 
of a Three Story Steel Frame Using System Identification", Report No. EERC 
78/24, Earth Engineering Research Center, University of California, 
Berkeley, 1978. 

13. Hays C., "Inelastic Materials Models in Earthquake Response", Proc. 3rd 
Engineering Mechanics Specialty Conference, September 17-19, 1979. 

14. fleck, J.L., "Determining Models of Structures From Earthquake Records", 
EERRL 78-01, Caltech, Pasadena, 1978. 

15. Ray D., et al, "Sensitivity Analysis of Hysteritic Dynamic Systems; 

Theory and Application", Comp. Methods in Appl. Mech. and Eng., Vol. 14, 
1978, pp. 179-208. 

16. Masri, S.F., and Caughey, T.K., "A Nonparametric Identification Technique 
for Nonlinear Dynamic Problems", ASME, Journal of Applied Mechanics, Vol. 
46, June 1979, pp. 433-447. 

17. Hudson, D.E., "Dynamic Tests of Full Scale Structures", Earthquake 
Engineering, Prentice Hall, pp. 127-149. 

18. Lee, Y.W., and Schetzen, J., "Measurement of the Weiner Kernels of a 
Nonlinear System by Cross Correlation", International Journal of Control, 
Vol. 2, 1965, pp. 237-254. 

19. Wiener, N., "Nonlinear Problems in Random Theory", M.I.T. Press, 

Cambridge, Mass., 1958. 

20. Udwadia, F.E. and Marmarelis, P.Z., "The Identification of Building 
Systems. II The Nonlinear Case", Bull, of the Seismological Soc. of 
America, Vol. 66, No. 1, February, 1976, pp. 125-151. 


37 



21. Udwadia, F.E. and Marmarelis, P.Z., "The Identification of Building 
Systems. II The Nonlinear Case", Bull, of the Seismological Soc. of 
America, Vol. 66, No. 1, February, 1976, pp. 153-171. 

22. Seed, H.B. and Idriss, I.M., "Analysis of Soil LiquiFaction: Niigata 

Earthquake", ASCE, SM3, May 1967, p. 83-108. 

23. Buzdugan, G.H., and Rades, M., "Forty-fourth Euromechanics Collequium on 
Dynamics of Machine Foundations", Bucharest, 29-31 Oct., 1973. 

24. Hadjian, A.H., "Support Motions for Mechanical Components During 
Earthquakes", Proceedings of the I. Mech. E. Conference, Nov. 7-9, 1978, 
pp. 27-46. 

25. Hitchings, D. , and Beresford, P.J., "A Comparison of Numerical Methods for 
Seismic Design of Mechanical Systems", I. Mech. E. Proceedings, November 
7-9, 1978, pp. 129-138. 

26. Kajimura, Y., and Shiraki, K., "Statistical Method of Estimating the 
Seismic Response of Light Secondary Systems", I. Mech. E. Proceedings, 
November 7-9, 1978, pp. 207-214. 

27. Hamming R., Numerical Methods for Scientists and Engineers, McGraw Hill. 


38 



TABLES AND FIGURES 


39 



Table 1. Description of linear system 
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Table! 2a. The coefficients of the identified stiffness of System 1 


i 

R, = b]y 

N/S 

R. = bj + bjy + + bjy^ + bjy'^ 

e . 

1 




Si 

Si 


s; 

1 

50 

0.001 

-0.00197 

48.130 

0.59691 

2.7956 

-1 .0551 

0.01552 

0.010 

-0.00513 

48.220 

0.35053 

2.5490 

-0.57827 

0.01527 


-0.01111 

48.208 



-0.19538 

0.01617 

2 

75 

0.001 

-0.01824 

74.250 

0.23820 

0.39391 

-0.13063 

0.00409 

0.010 

0.06529 

74.325 

-0.07721 

0.36444 

0.06970 

0.00442 

0.020 

0.17733 

74.335 

-0.44445 

0.36817 

0.27012 

0.00607 

3 

100 

0.001 

-0.00021 

99.436 

0.09831 

0.47977 

-0.10275 

0.00232 

0.010 

0.06144 

99.631 

0.01514 

0.14904 

0.11477 

0.00291 

0.020 

0.09254 

99.732 

0.03808 

-0.10082 

0.25814 

0.00452 

4 

50 

0.001 

-0.06270 

48.997 

0.26450 

0.62990 

-0.01724 

0.00906 

0.010 

-0.04072 

49.149 

0.48949 

0.51912 

-0.30196 

0.00896 

0.020 

-0.01943 

49.240 

0.68484 

0.44278 

-0.56636 

0.01198 
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Table 2b. The coefficients of the identified damping of System 1 


i 

HI 

M / Q 

D . = dj + d ] y + d 2 y ^ 


«i 

IN/ o 


^1 


-3 

1 

0.5 



0.51512 

- 0.00034 

0.00008 

0.04534 

0.010 

- 0.07826 

0.53241 

0.00085 

- 0.00006 

0.05553 

0.020 

- 0.07995 


0.00111 

1 

- 0.00012 

0.06109 

2 

0.4 

0.001 

0.03075 

0.41320 


0.00000 

0.03410 

0.010 


0.43327 

0.00008 

- 0.00002 

0.06786 




0.00005 

- 0.00003 

0.10025 

1 

0.3 

0.001 

- 0.00911 

0.30483 

- 0.00001 

0.00001 

0.02608 

0.010 

- 0.00942 

0.28528 

0.00057 

0.00009 

0.03947 

0.020 

0.02597 


■ 0.00102 

0.00016 

0.07600 

1 

1.0 

0.001 

- 0.15991 

0.95647 


0.00218 

0.03891 

0.010 

- 0.35716 

0.97538 

0.00655 

0.00146 

0.07359 

0.020 

- 0.50936 

0.98855 

0.00550 

0.00083 

0.11056 
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Table 3. Description of the system with nonlinear stiffness 

and linear damping 
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Table 4a. The coefficients of the identified stiffness of System 2 


1 

R .^ bjy + bjy ^ 

N/S 

R . 

1 

= bj + bjy + b^y2 + b^3 + t 






Si 

Si 

^3 


■ 



0.001 

0.00902 

48.915 


11.888 

- 0.42176 


1 

50 

10 

0.010 

0.02193 

49.039 

- 0.27482 


0.31841 

0.00809 

1 



0.020 

0.03476 

49.146 

- 0.76921 


1 .0329 

0.00815 




0.001 

- 0.02578 

74.441 

0.02301 

25.556 



2 

75 

25 

0.010 

0.06611 



25.566 

0.03328 

0.00365 




0.020 

0.18139 

74.490 


25.596 






0.001 

- 0.00499 

99.677 

- 0.02814 




3 

100 

20 

0.010 

0.10660 

100.160 

1.0550 

19.993 






0.020 

0.17982 

100.56 

2.2282 


- 3.0004 

0.00701 

■I 



0.001 

- 0.02028 

49.627 

0.02107 

- 9.7352 

0.08131 

0.00458 


50 

-10 

0.010 

- 0.02244 

49.146 

0.28733 


- 0.30165 

0.01273 

II 



0.020 

- 0.05730 



48.607 

0.60696 

- 9.1306 

- 0.71023 

0.02433 
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Table 4b. The coefficients of the identified damping of System 2 


i 

D , = bjy 

N/S 





^0 

"•l 

dg 

^3 

1 

0.5 

0.001 


0.52479 

- 0.00008 

0.00005 

0.061573 



0.52264 


0.00003 

0.05254 

0.020 

- 0.00323 

0.52365 

0.00115 

- 0.00002 

0.04916 

2 

0.4 

0.001 

- 0.01897 

0.40985 

- 0.00056 

0.000004 

0.04179 

0.010 

- 0.01952 

0.42118 

- 0.00070 

- 0.00001 

0.05485 

0.020 

0.01377 

0.42822 

- 0.00091 

- 0.00002 

0.06562 

3 

0.3 

0.001 

0.00748 

0.30273 

0.00001 

0.00003 

0.03412 

0.010 

- 0.02072 

0.30118 

0.00049 

0.00004 

0.03714 

0.020 

0.08008 

0.29381 

0.00043 

0.00006 

0.05204 

1 

1.0 

0.001 

0.08478 

1 .18900 

- 0.01297 

- 0.00428 


0.010 

- 0.30175 

1 .09520 

- 0.00124 

- 0.00217 


0.020 

- 0.64619 

1.00790 

0.00849 

- 0.00028 

0.10322 
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Table 5. Description of the system with both nonlinear stiffness and damping 


SYSTEM 3 (Stiffness = nonlinear) 
(Damping = nonlinear) 
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Table 6a. The coefficients of the identified stiffness of 
the second nonlinear system done by Method I and Method II 
in single precision and double precision calculation*. 


i 

Ri=bi] 

^+bsy^ 

N/S 

Method 

— 

Preci- 

sion 

R. = 
1 

. A , 

bo + bty 

A , /- 

+ bay + t 

'i 3 
>3y 


•H 

b3 

bo 

mm 

mm 

■B 


1 



■ 

I 

s 

-.01686 

48.375 

-.04161 

12.375 

.01429 

■ 




I 

D 


48.379 

-.04171 

12.370 


■ 




II 

S 


48.375 

-.04161 

12.375 


■ 



■ 


D 

-.01670 

48.379 

-.04171 

12.370 

.01426 

■ 




I 

S 

-.04644 

48.628 


11.606 

.01222 

B 

50 

10 

.01 

I 

D 

-.04626 

48.631 

-.02362 

11.601 

.01219 

■ 




II 

S 

-.04644 

48.628 


11.606 

.01222 





II 

D 

-.04626 

48.631 


11.601 

.01219 

B 




B 


-.07220 

48.700 

-.03593 

11.109 

.01391 

B 




H 


-.07202 

48.703 

-.03615 


.01388 

B 






-.07220 

48.700 

-.03593 

11.109 


B 





D 

-.07202 

48.703 

-.03615 


.01388 





I 

S 

-.03992 

73.254 

.10587 

22.916 

.00962 






D 

-.04018 

73.259 

.10634 

22.911 

.00960 






S 

-.10122 

71.991 

.22105 

25.168 

.01713 






D 

-.10089 

71.997 

.22135 

25.159 

.01710 






S 

.05966 

73.669 

-.06598 

22.185 

.00734 

2 

75 

20 



D 

.05936 

73.674 

-.06540 

22.180 







S 

-.05078 

72.433 

-.14566 

24.395 







D 

-.06778 

72.442 

-.14070 

24.390 

.01455 





I 

S 

.15867 

74.045 

-.21514 

21.440 

.00538 




.02 

I 

D 

.15835 

74.049 

-.21454 

21.435 

.00536 





II 

S 

.02881 

72.884 

-.50305 

23.505 

.01180 





II 

D 

.07893 

72.884 

-.51482 

23.486 

.01165 
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Table 6a. (Continued) 


1 



Method 

Preci- 

sion 

/N / 

R. = 

1 

■'i , C± 
30 + biy 

+ bay + bay 

E. 

1 

IB 

m 

N/S 

/\ . 
bo 

BB 

wmm 

mm 





H 

s 

-.01427 

98.843 

-.02418 

29.512 

.00446 





H 

D 

-.01415 

98.844 

-.02463 

29.513 

.00446 






S 

-.07047 

97.476 

.38302 

38.651 

.01348 






D 

-.07339 

97.477 

.38424 

38.656 

.01349 





n 

S 

.06645 

98.444 

,67074 

31.224 

.00772 

3 

100 

25 


I 

D 

.06663 

98.444 

.66983 

31.226 

.00772 





II 

S 

.03765 

97.113 

1.37590 

41.166 

.01813 






D 


97.110 

1.38350 

41.184 

.01787 





B 


.14580 

97.925 


33.115 

.01254 







.14605 

97,925 


33.120 

.01254 





H 


.10024 

96.607 

2.38390 

43.847 

,02338 





B 

D 

.19245 

96.610 


43.848 

.02465 

11111111 




B 

S 

-.03061 

49.290 

.07580 

-9.513 

.00789 




.001 

Bj 

D 

-.03054 

49.289 

.07600 

-9.512 

.00789 





B 

S 

-.05281 

48.922 

.23406 

-9.002 

.01204 

H 




B 

D 

-.05532 

48.922 

.23419 

-9.002 

.01203 

H 




I 


-.05748 


.05694 

-9.579 

.00656 


50 

-10 


I 


-.05664 


.05627 

-9.578 

.00657 

H 




II 


-.03826 

49.084 

.36383 

-9.207 

.01192 

H 




II 

D 

-.09496 

49.083 

.36194 


.01125 

H 




I 


-.08937 

49.501 

.04032 

-9.638 

.00586 

H 



.02 

I 


-.08861 

49.500 

.03962 

-9,638 

.00586 

|H 




II 


-.07520 

49.197 

.49163 

-9.367 

.01292 

1 




II 

D 

.10552 

49.199 

.49369 

-9.368 

.01675 


S - Single precision computations 
D = Double precision computations 
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* I = Method I 
II = Method II 

































Table 6b. The coef f icientij of the identified damping of the 

second nonlinear system done by Method I and Method II in 
single precision and double precision calculation*. 


1 

Di=di 

y+dsy^ 


Method 

Precision 

A A. 

D. = d 
1 

i 

0 + diy - 

A . A . 

1 • 1 • 

d2y + dsy 

3 


di 

4 

N/S 

A . 

do 

H 

A , 

d2 

d3 





I 

S 

--.00474 

.36882 

.00129 

.01888 

.04277 





I 

D 

-.00401 

.37062 

.00124 

.01881 

.04320 





II 

S 

-.00474 

.36882 

.00129 

.01888 

.04277 






0 

-.00401 

.37062 

.00124 

.01881 

.04320 





I 

S 

-.05902 

.35145 

.00412 

.02035 


1 

.3 


.01 

I 

D 

-.07636 

.34857 

.00502 

.02048 







S 

-.05902 

.35144 

.00412 

.02035 





. . .. 


D 

-.07636 

.34857 

.00516 

.02048 







S 

-.08631 

.33128 

.00487 

.02197 

.08588 






D 

-.03342 

.34376 

.00209 

.02147 

.08450 






S 

-.08632 

.33128 

.00487 

.02197 

.08588 






D 

-.03342 

.34376 

.00209 

.02147 

.08450 






S 

.04107 

.22486 

-.00077 

.02030 

.03131 







.03890 

.22599 

-.00073 

.02029 

.03146 






S 

.07877 

.21373 

-.00164 

.02056 







D 

.07516 

.21539 

-.00157 

.02053 

.03448 





I 

S 

.10625 

.23498 

-.00293 

.02029 

.03932 

2 

75 

20 


I 

D 


.23600 

-.00292 

.02028 

.03947 





II 



.23262 

-.00465 

.02052 

.04806 







.13476 

.23579 

-.00444 

.02048 

.04838 





I 

S 

.14124 

.26213 

-.00435 


.04924 




.02 

I 

D 

.13950 

.26321 

-.00431 

■■ 






II 

S 

.15053 

.27719 

-.00612 

.02004 






II 

D 

.19283 

.27513 

-.00666 

.02005 

.06244 
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Table 6b. (continued) 





N/S 

Method 

Precision 

II 

•H 

< n 

''i ■^i* '^i* '' 

do d- diy + d 2 y + c 

i* 3 
37 

*^1 

■ 


/\ , 
do 

■ 

d2 

A . 

dl 






S 

.04592 

.22146 


.02003 

.03488 






D 


.22193 


.02024 

.03505 






S 

.08132 

.26880 


.01952 

.06527 






D 

.08153 

.26875 

-.00566 

.01953 

.06539 





I 

S 

-.12147 

.16867 


mam 

.04661 

3 

.2 

.02 

.01 

I 

D 

-.11962 



mmm 

.04644 





II 

S 

.05635 

.24278 


.02032 

.06278 





II 

D 

.02908 

.24147 


.02034 

.06287 





H 


-.02790 

.13032 



.07451 




.02 


D 

-.27858 

.13092 

.00515 


.07438 






s 

.00825 

.22404 


.02094 

.06332 






D 

.10102 

.22733 

-.00345 

.02091 

.06714 





I 

S 

-.13261 

.22245 

.00044 

.01800 

.04685 




.001 

I 

D 

-.11423 

.22502 

-.00023 

.01790 

.04670 





II 

S 

-.09444 

.20412 

.00459 

.01831 

.04304 





II 

D 

-.07378 

.20734 

.00374 

.01822 

.04266 





I 

1 

s 

-.24066 

.22469 

.00251 

.01820 

.05645 

4 

.1 

.02 

.01 

I 

D 

-.22695 

.22739 

.00199 

.01813 

.05628 





II 

S 

-.00806 

.19753 

.00824 

.01835 

.05739 





II 

D 

-.04552 

.20146 

.00749 

.01825 

.05326 





I 

S 

-.28979 

.23868 

.00268 

.01811 

.06659 




.02 

I 

D 

-.28567 

.23928 

.00225 

.01809 

.06649 





II 

S 

.03103 

.20477 

.01147 

.01809 

.07794 





II 

D 

.21830 

.20312 

.01119 

.01810 

.09804 


*I = Method I S = single precision computations 

II = Method II D = double precision computations 
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Table 7. Comparison of the identified coefficients of the second nonlinear 

system in separable restoring force case by the separable restoring 
force (s) and the general restoring force (G) in Method I and single 
precision calculation. 



K=biy+biyH 

i • 5*9 

diy+day^ 


R.F. 


* 1 ^ _ 1 _ i J 

^9/ " ^0 + aiy -t 

i 2 , is, i* , i* 

Say + asy + + asy^ 

+ aeyj^' 

n. 

1 

■ 

i 

■ 


■ 

■ 

■ 

bi 

bf 

di 








S 

.02160 

48.375 

-.04160 

12.375 

.36882 

.00129 

.01888 

.04510 






.001 

G 

.00021 

49.997 

-.00728 

9.970 

.29566 

-.00014 

.02067 

.00024 

1 

50 

10 

.3 

.02 


S 

.10546 

48.628 

-.02344 

11.606 

.35145 

.00412 

.02035 

.06152 






.01 

G 

.01535 

49.937 

-.99801 

7.543 

.37732 

-.00208 

-.0043 

.05756 







S 

.15850 

48.700 

-.03590 

11.109 

.33128 

.00487 

.02197 

.08700 







G 

.07649 

50.552 

-5.9160 

-4.168 

.44665 

-.02794 

-.02717 

.11509 







S 

.00115 

73.254 

.10587 

22.916 

.22486 

-.00077 

.02030 

.03275 







G 

.00098 

74.999 

-.15306 

20.054 

.19541 

-.00064 

.02026 

.00011 

2 

75 

20 

.2 

.02 


S 

.16591 

73.669 

-.06598 

22.185 

.23498 

-.00292 

.02029 

.04000 






.01 

G 

.01562 

74.562 

-.29300 

21.541 

.21931 

-.00395 

.01926 

.03259 







S 

.29991 

74.045 

-.21514 

21.440 

.26213 

-.00435 

.02001 

.04953 






.02 

G 

.02943 

73.992 

-.42824 

22.887 

.24299 

-.00703 

.01801 

.06518 







s 

.03165 

98.843 

-.02418 

29.512 

.22146 

-.00128 

.02025 

.03517 






.001 

G 

.00035 

100.01 

.04555 

2i^.785 

.19621 

-.00036 

.02087 

.00005 

3 

100 

25 

.2 



S 

.05502 

98.444 


31.224 

.16867 

.00199 

.02174 

.04725 






.01 

G 

-.10294 

99.706 


24.804 

.29445 

.01986 

.01823 

.03422 







S 

.13320 

97.925 


33.115 

.13030 

.00516 

.02279 

.07556 







G 

-.19759 

99.228 

1.89125 

24.704 

.37880 

.03640 

.01822 

.06835 







S 

.16322 

49.290 

.07579 

-9.512 

.22245 

.00044 

.01797 

.04750 






.001 

G 

-.00004 

49.999 

.00324 

-10.00 

.09502 

-.00044 

.02033 

.00047 

4 

50 

-10 

.1 

.02 


S 

-.29814 

49.405 

.05694 

-9.579 

.22469 

.00251 

.01820 

.05683 






.01 

G 

.01860 

50.085 

.16231 

-10.352 

-.00928 

-.03925 

.01660 

.03874 







S 

.37916 

49.501 

.04032 

-9.638 

.23868 

.00240 

.01811 

.06685 






.02 

G 

.03456 

50.156 

.33714 

-10.732 

-.10973 



-.07193 

.11261 

.07734 
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Table 8. The exact parameters and the identified corresponding terms 
used in the general restoring force case. 


£ 

N/S 


= ^£^£ + ^£^£ + ^ 

a/TvTO A/A/ A/A/ A/A/A/ 

n 


^£ 

^£ 


®£ 

^£ 


EXACT 

50. 


10. 

.3 

.02 

0.0 



.001 

49.961 

9.743 

10.010 

.29432 

.01974 

.105 

.00423 

1 

.01 

50.083 

1.067 

10.014 

.17897 

.01090 

15.467 

.04473 


.02 

48.998 

5.400 

9.9320 


.01350 

33.825 

.09205 


EXACT 

75. 

20. 

0.0 

.2 

.02 

.2 



.001 

75.156 

17.030 

-.0617 

.1925 

.0211 

2.8926 


2 

.01 

75.493 

17.097 

-.0819 

.2172 

.0168 

1.6016 



.02 

74.262 

49.431 

-.1786 

.2426 

.01302 

.7544 

.03081 


EXACT 

100 

25 

10.0 

.2 

.02 

0.0 



.001 

99.946 

21.343 

9.9850 

.32668 

.00575 

o 

o^ 

.00094 

3 

.01 

97.829 

44.649 

10.671 

.37731 

.01309 

3.53 



.02 

88.006 

25.434 

-2.2157 

.47349 

.01451 

11.718 



EXACT 

50 

-10 

0.0 

.1 

.02 

.2 



.001 

49.949 

-9.990 

.03473 



.19440 

.00136 

4 

.01 

49.484 

-9.870 

.35721 


.10968 

.13410 

.01341 



.02 

48.997 

-9.879 

.6953 

-.21806 

.19722 


.02677 
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Figure 1. Close-coupled system 
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Figure 3. Four-degree-of-f reedotn example 
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Figure 7. A comparison of the response of the actual system ( — ) 

with that of the identified system ( — — — ), (a) under 
base excitation, (b) under an impulse force of the ten 
units applied at M for System 1. 
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Figure 9. A comparison of the response of the actual system ( — ) 

with that of the identified system ( — — — ), (a) under 
base excitation, (b) under an impulse force of the ten 
units applied at for System 2. 
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DISPLACEMENT VELOCITY 

Figure 12. A comparison of actual stiffness and damping 

denoted by solid ( — ) lines and the identified 
results by Method I and single precision ( — — ), 

Method I and double precision ( ), Method II and 

single precision ( — - — ) and Method II and 
double precision ( — - - — ) in nonlinear case (l-bl 
of N/S = 0.00 
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DISPLACEMENT 


Figvire 13. A comparison of actual stiffness and damping 

denoted by solid ( — ) lines and the identified 
results by Method I and single precision ( — — ), 

Method I and double precision ( ), Method II and 

single precision ( — - — ) and Method II and 
double precision ( — - - — ) in nonlinear case (1-b) 
of N/S = 0.01. 













DISPLACEMENT 


A comparison of actual stiffness and damping 
denoted by solid ( — ) lines and the identified 
results by Method I and single precision ( — — ), 

Method 1 and double precision (-■ ), Method II and 

single precision ( — - — ) and Method II and 
double precision ( — - - — ) in nonlinear case (1-b) 
of N/S = 0.02. 
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Figure 16. A comparison of actual restoring force denoted 

by solid ( — ) line and the results identified by 

the general restoring force method ( ) for 

the general restoring force case with N/S = 0.01. 
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